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The heat kernel in curved space-time is computed to fourth order in a strict expansion in the 
number of covariant derivatives. The computation is made for arbitrary non abelian gauge and 
scalar fields and for the Riemann connection in the coordinate sector. The expressions obtained 
hold for arbitrary tensor representations of the matter field. Complete results are presented for the 
diagonal matrix elements and for the trace of the heat kernel operator. In addition, Chan's formula 
is extended to curved space-time. As a byproduct, the bosonic effective action is also obtained to 
fourth order. 
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I. INTRODUCTION 

The heat kernel operator, the exponential of the Klein- 
Gordon operator, is a useful tool in quantum field theory 
since it is ultraviolet finite, one- valued and gauge covari- 
ant and allows to obtain the propagator and the effective 
action P, [3, H, Quite different applications of the 
heat kernel (such as spectral densities, index theorems, 
C-function, quantum anomalies, chiral gauge theories, ef- 
fective theories of QCD, Casimir effect, black hole en- 
tropies, membranes, etc.) are illustrated in 0, H, 

i, El, El m El- The exact evaluation of the heat 
kernel is not possible in general (see however, [l^ [l^ 
and [13, [H, [3 particular cases). Nevertheless the 
series expansion in powers of the proper time is available 
with computable coefficients, the so-called Hadamard- 
Minakshisundaram-DeWitt-Seeley (HMDS) or heat ker- 
nel coefficients [1, [2^ . These coefficients have been com- 
puted with different techniques in several setups and to 
dimension ten [U, [13, HI, Hj, HI HI, [13 . For reviews on 
spectral geometry and field theory on curved space see 
[28, 29, 30, 31, 32, 33, 34, 35, 36]. For an alternative ap- 
proach based on covariant perturbation theory see 37 1. 
Extensions to finite temperature can be found in [3i,[3|. 



Exp ansions around non c-number mass terms are given in 
[ia,[ill[i3l- The extension to non commutative quantum 
field theory has been presented in [isj . 

In [44,] we investigated a different expansion for the di- 
agonal matrix elements of the heat kernel operator where 
the terms are classified by their number of covariant 
derivatives. This can be regarded as a resummation of 
the HMDS coefficients to all orders in the (non abelian) 
mass term, which needs not be small. There, explicit re- 
sults were presented for boundaryless compact flat mani- 
folds to four derivatives for the diagonal matrix elements 
and to six derivatives for the trace of the heat kernel 
operator, for Klein-Gordon operators with arbitrary non 
abelian gauge and scalar external fields. In the present 
work we extend those results to the case of curved man- 
ifolds endowed with the Levi-Civita connection in the 
coordinate sector. 



The scope and ideas involved as well as some defini- 
tions are given in Section [Hi In particular, we avoid the 
standard approach of reducing tensor fields (e.g., the pho- 
ton field) to scalars by introduction of a tetrad [sj . That 
approach allows to use the heat kernel formulas derived 
for coordinate scalars but at the price of complicating 
the internal space (due to the new tetrad index) with its 
corresponding modified connection. Instead, we provide 
formulas which are equally simple regardless of the tensor 
representation of the fields (on which the Klein-Gordon 
operator is acting) using only the original gauge connec- 
tion in the sector of internal indices and the Levi-Civita 
connection in the coordinate sector. 



In Section IIIII we discuss the shortcomings of Chan's 
approach in the curved case and develop our own ap- 
proach based on the use of covariant symbols [i^, [40 ]. 
The covariant symbols arc in fact multiplicative opera- 
tors (with respect to x) which define a faithful repre- 
sentation of the algebra of pseudodifferential operators. 
This technique is similar to that of symbols for pseudod- 
ifferential operators [H [13, [H, [4^ and shares with it the 
feature of providing diagonal matrix elements of generic 
operators, not necessarily related to the heat kernel, since 
it is not based on recurrence relations. However there is 
an important difference between both techniques: in the 
standard method of symbols covariance is not manifest 
prior to momentum integration whereas with covariant 
symbols covariance (under both gauge and coordinate 
transformations) is manifest at every step of the calcu- 
lation. In that section results are provided for the diag- 
onal matrix elements to four derivatives. These results 
are presented in compact form using the technique of la- 
beled operators, also used in |4j]. In this notation the 



mass term carries a label indicating its position in the 
expression; in this way it becomes effectively a c-number 
and momentum integrals can be carried out explicitly. 

Explicit formulas for the derivative expansion of the 
trace of the heat kernel operator are presented in Section 
IIVI This is of interest in the computation of the effective 
action. It it is noted that many, in fact most, of the new 
terms introduced by the curvature can be eliminated by 
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a suitable redefinition of the mass term in the Klein- 
Gordon operator. 

In Section |V] we obtain the expressions for the diago- 
nal matrix elements and trace of the heat kernel in the 
form first derived by Chan [5flj for the effective action to 
four derivatives and extended in [sij to six derivatives. 
This is done following a rather indirect path since Chan's 
derivation is not used. Instead a Chan's form is proposed 
and the coefficients are adjusted to reproduce the results 
previously derived using covariant symbols. The ques- 
tion of whether Chan's elegant method can be used in 
curved manifolds is left open. In any case, as expected 
from the experience in the flat case, Chan's form is truly 
much more compact. Drawbacks are that the expres- 
sions are less explicit because one parametric integral is 
left undone, and for the same reason there is a by parts 
integration ambiguity in the formulas. 

Section IVII is devoted to computing the bosonic effec- 
tive action in an explicit way, with the help of labeled 
operators. The necessary momentum integrals can be 
done with dimensional regularization. They are given in 
[S^ l using the MS scheme. 

Somewhat more explicit details on the calculation with 
covariant symbols are given in Appendix |^ Trace cyclic 
property and integration by parts identities are derived 
in Appendix [Bl 



internal space. Our convention for the Riemann tensor 
is such that, if t/i is a world vector. 



where 



[V^, V^ItZ-A = f^M>>V^A + Rf^^Xa^cT (2.1) 



^fiu = df^uj^ - di^LUf^ + [ujf,, Lu^] (2.2) 



is the gauge field strength tensor. In addition, for the 
Ricci tensor and scalar curvature 



(2.3) 



In the previous formulas we use the same notation for 
contravariant and covariant indices. We will follow this 
convention throughout unless an ambiguity arises. 
The Klein-Gordon operator is of the form 



K ^ .g^'V^V, + X . 



(2.4) 



X{x) is a scalar field with respect to world indices and 
a matrix with respect the internal space. As is well- 
known the heat kernel operator e'^^ is ultraviolet finite 
for Re (r) > 0, and its matrix elements admit an asymp- 
totic expansion in powers of t. For the diagonal matrix 
elements 



II. DERIVATIVE EXPANSION OF THE HEAT 
KERNEL 



We assume a compact boundaryless Riemannian man- 
ifold of dimension d and Euclidean metric g^^{x) to rep- 
resent the space-time upon Wick rotation. Further, for 
the world sector indices^ we take the Levi-Civita connec- 
tion (i.e., torsionless and metric preserving). The Klein- 
Gordon operator will act on wavefunctions il>{x) (matter 
fields) defined on the space-time manifold. The wave- 
functions are vectors with respect to some representation 
of a certain gauge group. Without loss of generality we 
will assume that ip carries a single gauge (or internal) 
index, ij) is also allowed to carry world indices, that is, 
we do not assume V' to be a world scalar.^ Following 
for instance [34], we use a single covariant derivative 
which acts on all indices with the appropriate connection, 
V = d + V + u). is the connection on world indices 
and a;^(x) the connection on gauge indices, a matrix in 



1 



(47rT)'^/2 ^„ 



T"a„(a;) . 



(2.5) 



The coefficients a„ are still operators with respect to 
gauge and world indices since the brackets (a;|, \x) re- 
fer only to x-space. To lowest orders, the well-known 
result is 



ao = 1 
ai = 

1 



X + -R, 
6 



^2,1^ 1 ^2 1^,. 1 



2 6 12 ^'^ 6 30 '^'^ 



72 180^ 



180 



-R 



(2.6) 



In these formulas (and hereafter) we use the following no- 
tational convention: the covariant derivative of an object 
is represented by adding a world index to it and further 
covariant derivatives add further indices to the left. So 
for instance. 



In this work we will use the label world interchangeably with co- 
ordinate or space-time in expressions like "world tensor" , "world 
index", etc, to refer to properties tied to indices fi,u,..., as- 
sociated to natural bases, d/dx'^ , of the tangent space of the 
Riemannian manifold. 

Of course one can choose to transform such world indices into 
internal or gauge indices using a tetrad field. This is the standard 
approach [SJ]. Our formulas hold whether this choice is made or 
not. 



A^ = [V^,A], A^, = [V^,[V.,A]], 

We have also introduced the quantity 

:= [V^,V,]. (2.8) 

Let us emphasize that Zfj,^, is different from r?^,^. The op- 
erators Zfj^iy and $1^1/ are both multiplicative with respect 
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to a;-space and matrices in internal space, however, un- 
like ri^jy, Z^i, acts also on world indices (cf. (|2.ip ). This 
implies, for instance, that Z^^, and T^a^ do not commute: 

\Z^v, TZap] = Ftfj.i'aaTZap + Rfiv/SaT^aa ■ (2-9) 



Likewise, [Zfj,^,Xx] = [Qf,i,,Xx] + R^j^^xaXa, etc. The 
derivative of Z^j^^, will also be needed,'^ 



dimension, 
A2 = 



1 



6 



-RX 



12 ^"^^ 30 '^'^ 

I \_ 7->2 I 

180 '^""'^ 



72 



— 7^2 

180 



(2.12) 



(2.10) 



In our formulas we will use Z^^, rather than fJ^j/. Nev- 
ertheless, if desired one can always move the Z's to the 
right using their commutation properties and apply them 
to the wavefunction ijj to produce fi's and i?'s. (Of 
course, the result so obtained will depend on the concrete 
tensor representation of the wavefunction.) For instance, 
the term Z^^/12 in a2 is equivalent to more standard 
when the wavefimction happens to be a world 
scalar. The advantage of using Z^^...^^ is that the expres- 
sions take the same form regardless of the world-tensor 
representation of ^ and yet the original connections T 
and uj are used. (See footnote [H) Fuller details can be 
found in [46j . 

By dimensional counting it is clear that the standard 
heat kernel expansion (|2.5p is an expansion with coef- 
ficients ordered by their mass dimension: if g^i^ carries 
dimension 0, carries dimension 1, and X carries di- 
mension 2, the coefficient a„ has dimension 2n. In what 
follows we will set t — I since dimensional counting al- 
lows to restore the parameter r at any time if needed. 

In this work we consider a different classification of 
terms in the heat kernel, namely, operators are classified 
by the number of covariant derivatives they carry, so 
counts as order 1 while X and g^^ count as order cero. 
Thus, e.g., XX/^t is of first order, R^va/i and of second 
order, Zani, of third order, and so on. In this expansion 



{x\e 



oo 



(47r)<i/2 

^ ' rt=0 



An(x) , 



(2.11) 



where the coefficient A„ collects all operators with 2n 
derivatives and any number of X. In turn, the An can 
be reexpanded using operators classified by their mass 



^ This is an exception to the index convention. The second term 
in the definition of Zapuv is required to make it a multiplicative 
operator with respect to x. E.g. 

Similar terms appear in higher derivatives, Za^pfiu, etc, [4^ . 



(Of course, when r is not set to unity these coefficients 
are also functions of r.) 

In [iJl we presented complete formulas (i.e., valid to 
all orders in X) for j4o. Ax and A2, as well as Sq, -Si, 
and ^3 (defined below), for the flat case. Presently we 
extend those results (except B^) to curved manifolds. 

The key tool to write the result to all orders in X in 
closed form is the use of labeled operators. For instance, 
all terms in Ai of the type X'^X^^X"^^ with n,m = 
0, 1, . . ., can be collected as 



I2,2X^^ 



(^1 



(^1 



X, 



(2.13) 



The label 1 in Xi indicates that the corresponding X 
should be placed just before (i.e. to the left of) the fixed 
operator Xf^f^, likewise, X2 indicates that that X is to 
be put just after (i.e., to the right of) the fixed operator. 
Upon a series expansion in powers of Xi and X2 



l2,2Xnn 



1 

12 



1 

12' 



:X2 



+—X^ + —X: 

An ^ An 



40 



40 



30 



X1X2 



X„ 



TXnn + —XXun + —XnuX 



12 



12 



1 2 1 2 

40 40 



30 ^'^ 



(2.14) 



Likewise, X^^XaXaX^ would stand for X'^X^.XX^X, etc. 
The labels always refer to the position of X^s with respect 
to "fixed operators" such as Xafj.^, Zf^,^, etc. Because X 
commutes with R^i,ap we will not need to include the 
Riemann tensor among the set of fixed operators. So for 
instance, l2,2RX^^ will used to mean R multiplied by 

l2,2X^^. 

The important point is that the labeled operators can 
be treated as c-numbers, e.g., X1X2 — X2X1; the true 
position of the labeled operator is given by its label. This 
is similar to what happens within normal or chronological 
orders. Labeled operators were introduced in [SSl [53|. 

The function I2.2 belongs to the family of functions 
Iri,...,r„, with arguments Xi, . . . They are defined 

as 



dz 



(2.15) 
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where 



N = {z-xy 



(2.16) 



and Ni = {z — Xi)^^ , and F is a positively oriented closed 
simple path on the complex plane enclosing the eigenval- 
ues of X (at the given point x). These functions enjoy 
a number of properties described in ^4^. In particular 
they are entire functions of the Xi and can be computed 
by recurrence relations starting from Ii = e^'^ . They are 
not linearly independent; integration by parts implies the 
relation 



-^ri,r2,. ..,!•„ — '^^Tjlrx,r2,...,rj + l,...,r 



The trace of the heat kernel operator, 



Tr(e-^) = j d'^x^ty:{x\e'^\ 



(2.17) 



(2.18) 



is also of great interest in applications, such as the com- 
putation of the effective action (see Section IVip . The 
trace tr refers to gauge and world indices. The trace 
of the heat kernel admits an asymptotic expansion with 
terms classified by their mass dimension: 

1 °° r 

^ JI^JlJ d'^Vg^^Kix)), (2.19) 



the coefficients a„(a;) are local, i.e., they depend on a fi- 
nite number of derivatives of the external fields (namely, 
the metric, the gauge connection and the scalar field X). 
Consequently, the expansion at a given point xq is not 
sensitive to modifications of these external fields taking 
place outside a fixed neighborhood of xq. Such modifi- 
cations affect the exact matrix element at xq (and hence 
its T dependence) but this is not seen by the asymptotic 
expansion. By the same token, the derivative expansion 
is expected to be asymptotic too, since the coefficients 
An{x;T) are also local. To make mathematical sense of 
the asymptotic series it is necessary to write it as a power 
series expansion. The derivative expansion can be viewed 
as a power series as follows: for a given point xq consider 
a deformation of the external fields such that within a 
fixed neighborhood of xq the deformation is just a di- 
latation by a parameter A 

X{x)^X{xx), 
uj^{x) 1-^ Xujf_,{xx), 



= + X{x^ - x^) 



(2.23) 



The deformation is smoothly continued outside of the 
neighborhood. This produces a family of Klein-Gordon 
operators -?Ca-* By construction the parameter A counts 
the number of derivatives, that is, 



with 



bo 
hi 



72 



1 



1 



180 



180 



Likewise, the terms can be classified by their number 
of derivatives 



Tr(e^) 



1 



^ ' n=0 



d''xy^tT{B„{x)). (2.21) 



Due to the trace cyclic property and integration by parts, 
there is an ambiguity in the definition of Bn{x). This 
allows to bring Bn{x) to a simpler form starting from 
An{x). In turn, functional derivation allows one to obtain 
An from Br, 



An{x) = 



5X{x) 



d'^Xy^tT{Bn{x)). 



(2.22) 



Before closing this Section, let us comment on the na- 
ture of the derivative expansion in the present context. 
[We reinstate r for this discussion.] As is well-known, 
the standard expansion in powers of r is an asymptotic 
one, reliable for small r only. Intuitively, this is because 



Anixo;T,X) = A^"A„(a;o;r) 



(2.24) 



Therefore the derivative expansion can be read off from 
an expansion of {xo\e'^^^\xo). This construction suggests 
that, although the series in A is asymptotic, the coeffi- 
cients An{x; t) themselves are well defined quantities, as 
is also the case for the coefficients an{x). This point is 
to be settled by a rigorous mathematical approach. Note 
that this definition of the coefficients (from an expansion 
in A) is not identical to the alternative definition 



Anix;T) 



[xV 



(2.25) 



(using am,n{x) to denote the terms of am{x) with 2n 
derivatives). We expect both definitions to coincide. 
This expectation relies on (i) the fact that the functions 
Iri,r2,...,r„ are entire functions of Xi and so of r, thus 
their expansion in r is absolutely convergent, and (ii) 
the obvious reason for the series in r to be asymptotic 
does not apply here since An{x; r) is itself local. In sum- 
mary, we expect the derivative expansion to be valid for 
finite (not small) r, provided A is sufficiently small. 



* The dilatation depends not only on xo but also on the coordinate 
system used. Eq. 112.241 1 is coordinate independent. 
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III. DIAGONAL COEFFICIENTS 

In [i^l the calculation of An in flat space-time was 
based on that of i?„. In turn i?„ was adapted from the 
result of Chan [50] for the effective action to four deriva- 
tives and its extension in [5l| to six derivatives. Unfortu- 
nately, it is not obvious how the elegant approach of |50| 
is to be extended to the curved case. As we may recall, 
in that approach a symbol method (Z?^ — > 13^ + Pp.) is 
applied, to Tr(log(_ftr)). The expression is then formally 
expanded in the number of derivatives and brought to 
a canonical form by using integration by parts (with re- 
spect to Pfj_) as well as formal cyclic property. This canon- 
ical form is not manifestly gauge invariant but it allows 
to unambiguously identify the originating gauge invariant 
expression from which it comes, because no two different 
gauge invariant expressions may have the same canonical 
form. The trouble is that a similar statement does not 
hold in the curved case. For instance. 



(3.1) 



(with N — {p^ — X)^^), since p^ and X commute. How- 
ever, under formal cyclic property the same expression 
would be equivalent to 

Tt{p^N'^[Z^„,P^]) = Trip^N^R^^^^p^) 

= Tiip^p.n^^N^) . (3.2) 

This is equivalent to —TT{p^RN'^)/d and does not vanish. 

Since Chan's approach is not available, we will com- 
pute An from scratch and then obtain Bn from it. The 
starting point is the representation 

(3.3) 



p 27ri z 



X 



where we can apply the method of covariant symbols |45l . 
l46j . This is the second key ingredient of our approach. 
This gives 



{x\e 



dz 1 
r 27ri ^ 



d'^p 



— ■ (3.4) 
X 



The covariant symbols are designed for computing diag- 
onal matrix elements of general operators, no necessarily 
the heat kernel. The covariant symbols were introduced 
by Pletnev and Banin in [4^ for the gauge connection 
and extended to curved space-time in [46[ where they 
have been computed to four derivatives. Explicitly, to 
two derivatives 



X 

=2 



X-X^d" 
o 1 



1 

2!' 



XalS dad. 



'/3 



1 



Pj, + g-R + ZaliPadp - -TZaf3 Padf3 
1 



RafiPuPaP/sd^idi, 



(3.5) 



where 9^ = d/dp^ and the dots refer to terms with three 
derivatives or more. Note that, unlike ordinary sym- 
bols, these operators are multiplicative and covariant al- 
ready before momentum integration. The calculation is 
straightforward along the lines of the examples presented 
in [46]. Details are provided in Appendix One obtains 



Ao = Ii=e 



X 



(3.6) 



to zeroth order and 



Ai = (/i,2 - 2/i,3)X^p + (2/i.i,2 - 4/i,i,3 - 2/1,2,2)^, 



-\hR 



(3.7) 



to second order. Using the identities (|2.17p . this can be 
rewritten as 



Al — l2,2Xaa + 2/2 1 2^^^ 



6 



(3.8) 



Let us emphasize that the functions /(A"i, ^2) multiply- 
ing and f{Xi,X2,X:i) muhiplying X^, are well de- 
fined and unambiguous. The only ambiguity enters in 
how they are written in terms of the overcomplete basis 

The two terms without R are identical to those of the 
flat case in [i^l, only with the "minimal coupling" re- 
placement D — d + uj^W = d + T + uj. These are 
then minimal terms required by covariance (under gen- 
eral coordinate transformations). The term with R is 
non minimal; covariant but not required by covariance. 
In general we will obtain coefficients of the form 



Ar^ 



a: 



At 



(3.9) 



The terms in An are those which contain explicitly the 
Riemann tensor and so vanish in the flat space case. On 
the other hand, the minimal terms A!^ can be recon- 
structed by minimal coupling from the flat space expres- 
sions. Let us warn, however, that this separation is not 
an unambiguous one: in general A™ will depend of the 
concrete An of flat space used to apply the minimal cou- 
pling. This is because reordering of world indices in the 
flat space expression may introduce Riemann tensors in 
the curved case.^ 

The calculation of A2 gives 
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For instance, 

I 



— 2/24,2 Z/^ii^Z^^ 

+ (2/2,2,2-4/3,0.3 + 4/2,1,3) 

+(2/2,2,2 -4/3,0,3 + 4/3,1,2) 

+ (4/2,2,1,2 " 16/3,0,1,3 — 8/3,0,2,2 + 8/3,1,1,2) X^Zfj^^Xjj 
+(16/3,0,1,3 — 8/2,1,1,3 — 4/2,1,2,2 + 8/3,0,2,2) Zfj_i,X^Xi, 
+ (16/3,1,0,3 — 8/3,1,1,2 — 4/2,2,1,2 + 8/2,2,0,3) Xfj^X^Zfj^i, 
+2/3,3 Xf^f^ 

+ (2/2,3,2 + 4/3,1.3 + 2/2,2,3 + 2/3,2,2) X^,^,X,, 

+ 8/3,1,3 X^yX^iJ 

+ (8/3,1,3+ 4/3,2,2) ^M'^.^M 
+ (8/3,1,3+ 4/2,2,3) 

+ (4/2,2,2,2 + 16/3,1,1,3 + 8/2,2,1,3 + 8/3.1,2,2) X^X^ijXy 
+ (2/2,2,2,2 + 8/3,1,1,3 + 4/2,2,1,3 + 4/3.1,2,2) Xf^iXy^Xf^i 

+ (2/2,2,2,2 + 8/3,1,1,3 + 4/2,2,1,3 + 4/2,3,1,2 + 4/3,1,2,2 + 4/3,2,1,2) X^^Xi,Xi, 
+ (2/2,2,2,2 + 8/3,1,1,3 + 4/2,2,1,3 + 4/2,1,3,2 + 4/3,1,2,2 + 4/2,1,2,3) X^X^Xi,i, 
+ (16/3,1,1,3 + 8/3,1.2,2) X^ijX^Xi, 
+ (16/3,1.1,3 + 8/2,2,1,3) X^X^Xi,f^ 

+ (4/2,2,1,2,2 + 16/3,1.1,1,3 + 8/2,1,2,1,3 + 8/3.1,2.1,2 + 4/2.2,2.1,2 + 4/2.1,2.2,2 

+ 8/2,1.3,1,2 + 8/2,2,1,1,3 + 8/3.1,1.2,2) X^X^XyXi, 
+ (4/2,2.1,2.2 + 16/3,1.1,1.3 + 8/2.2,1.1,3 + 8/3.1,1.2,2) Xf^iX^Xf^iXy 

+ (4/2,2.1,2.2 + 16/3,1.1,1.3 + 8/2.2,1.1,3 + 8/3.1,1.2,2) Xf^iX^X^X^ . (3.10) 
^2 — -^hR^l.^J^ + 72^^"^^ ^ Q^^:2RXf_,f_, + -l2,2RtiX^ + -h^zR^^^X^y + -/2,1,2-R-X'^ 

+ (^-^2,2,2 + 2^^2,3,2 - -^^3,1,3)'^A';'^M^"^ ^ 180^^^^" ^ 180^^^'^""'' ' (^-H) 



is formally identical to the expression in (4.10) of 
[4^ . (There mirror symmetry of An was exploited to 
write A2 in a shorter, less explicit form.) As said before, 
in the formula for A2 the Riemann tensor is not to be 
considered as one of the fixed operators. 

The coefficients Ai and A2, for curved space and non 
abelian gauge group, with contributions to all orders in 
X are computed here for the first time. Upon expansion 
in powers of X they reproduce the corresponding terms 
with up to four derivatives of the standard expansion a„. 
In particular the coefficients quoted in Eqs. (4.26-29) of 
[34] are correctly reproduced. (The comparison with [33 | 
is achieved by restricting our results to the case of world 
scalar wavefunction and using the trace cyclic property, 
but not integration by parts.) Results to all orders in 
X and up to four derivatives for curved space are also 
available from [s^ for the so-called minimal case, i.e., 



X{x) a c- number and tj^ = 0, and "tp a world scalar. Our 
formulas also check Eq. (22) of [ssj . 



IV. COEFFICIENTS FOR THE TRACE 



Application of integration by parts and the trace cyclic 
property allows to obtain the simpler coefficients /3„ . To 



7 



four derivatives they take the form 



with 



Bn = h 



6 



B2 — { — 12,2,2,2 + 'il3,l,3,l)XtJ.XfiX^Xi, 

+ -;^I2,2,2,2X^X„X^X^ + AI^ 1,^X^X1^1X^1, 

+h,zXfj_nXi,y + 2l2,2,2Xfj_XiyZfj_^ + -l2,2Z^ 

30 w 72 ^ 



1 

180 



(4.1) 



[For short we have written 12,2 for 12,2,0 , etc.] The sim- 
plest way to obtain this result is starting from its Chan's 
form, to be discussed in the next section. The vaUdity of 
integration by parts for general world representations of 
the wavefunction space is shown in the appendix. The 
cyclic property for multiplicative operators is also subtle 
when Z^i, (or more generally Z^^...^^) is present because 
this operator acts on world indices. For instance, under 
trace, 



X,/XpZpi, 



(4.2) 



The last term would not be present under the standard 
trace cyclic property of matrices. The rationale of the 
extra term is that in X^Z^^X^, Z^^^ acts on the index v 
of X^ (as well as on the subsequent world indices in the 
wavefunction) while in X^jX^Z^j^i, it does not; the missing 
contribution is added by the extra term. This equation 
is derived in Appendix [Bl^ 

It can be verified that An and _B„ are indeed equivalent 
inside /d''a;^tr( ) and that ([2:2^ is fulfilled. 

To this order one can see that most of the non minimal 
terms can be generated by a modified "minimal coupling" 
prescription. Namely, d + u ^ d + u) + T and X — > X', 
with the new scalar field 

X' = X + ^R+^{Rp,-nl,+Rl,^p)+0{V^) . (4.3) 

The trace of the heat kernel can then be expanded as 



1 °° r 



(4.4) 



Of course, relations like 1)4. 2p imply that starting from the same 
Bn of flat space but written in different ways will in general yield 
different results for the minimal coupling part of B„ in curved 
space. 



B[ 



2^2,2^1^1^ pi ' 



B2 - (-^2,2,2,2 + 4/34,3,1)^^^^^;!,^;!, 

+ 2-^2,2,2,2-'^^-'^^-'^^-'^^ + ^Iz,l,Z^'pi^'p.^'vv 
'^^L'i^'pilJ.^lv + 2^2, 2, 2-'^^-'^^-^^'' + 2^2,2ZpivZpiv 



(4.5) 



{^ri,...,r„ being defined as in (|2.15p but using X' instead 
of X.) There are corresponding coefficients A'^- The 
relations J d'^Xy^tT{An) = J dfix^U{Bn) and (P?^ 
hold also for the primed coefficients. The redefinition 
X' = X + ii? is quite standard in the literature [l5| to 
eliminate some of the terms. Note that beyond that the 
primed expansion is no longer a strict expansion in the 
number of covariant derivatives (and X' will depend on 
r when t is restored). 



V. CHAN'S FORM OF THE COEFFICIENTS 

Let us call the coefficients A„ and B^ just derived the 
coefficients in X-form, to distinguish them from their 
Chan's or A^-form, to be discussed in this section. 

Let us briefly summarize Chan's method in flat space 
jsoL IsH ] . In this method the results are obtained in terms 
of derivatives of7V= (z — A)~^, instead of derivatives 
of X, that is, Np = [V^A^], N^^, etc. Because z appears 
inside Nf^^...^j_^ the integral over z cannot be carried out 
explicitly and is left undone. Integration by parts (with 
respect to z) allows to reorder terms so that in each 
term of Bn the quantity N (derivated or not) appears 
exactly 2n times. A virtue of this approach is that for 
each Bn there is only a limited number of available co- 
variant structures constructed with 2n A^'s and 2n V's, 
thus the expressions so obtained are quite compact. To 
pass a result given in A^-form to X-form is, of course, 
straightforward using the relation A"^ = NX^N and its 
derivatives. As pointed out before, the result in A-form 
is free from ambiguities. 

Because the extension of Chan's method to curved 
space is not known, the existence of a Chan's form for the 
heat kernel coefficients in the curved case is not obvious. 
In principle, undoing the z integrals in the /ri,...,r„ and 
using identities of the type 



X„, = 



Xp, = N^^Np^N-^ 



(5.1) 



would allow to bring the coefficients computed in A"-form 
to an almost Chan's form, except that, in general, neg- 
ative powers of A^ will be present. Moreover the precise 
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result will be subject to ambiguities due to integration 
by parts on z and not algorithm is available to bring it 
to a compact form. This method works for and B\. 

The procedure that we have followed to obtain an A^- 
form for the coefficient Bi is as follows. In the flat case 
the A^-form is known for B^ and hence for A2 (by func- 
tional variation with respect to X), therefore we apply 
minimal coupling there. This already reproduces most of 
the terms of the known full Ai in X-form, (|3. 10113. lT|) . 
For the few remaining terms of A2 it is relatively easy to 
bring them to a rather compact A^-form by hand. This 
remainder is the functional variation of the non minimal 
remainder in i?2, i?2^, not yet determined. To obtain 
B^ we simply write down the most general terms hav- 
ing four derivatives, at least one Riemann tensor and no 
more than four A'^'s and with arbitrary numerical coef- 
ficients. These coefficients are then chosen to reproduce 
the non minimal remainder of A2. This procedure gives^ 



Bo 
B2 



, 1 



2 ^ 



(5.2) 



-2NN^N,NZ^, 



12^^^^ 



72 180 180 '""^^ 



where we use the shorthand notation 



dz 



(5.3) 



B2 in X-form, (14.11). has been obtained from this A^- 



form. Eqs. (15. 2|) are the extension of Chan's formulas to 
curved space-time. Once again, to this order, the only 
non minimal term surviving is —'^TZ^^NNf^NNi, if X' is 
used throughout. 



This gives® 

+2NNIn^^N + 2NN^^N^N 

+2N^Nf^^N + 2NNf,f,NN^^N + 2NN'^^N^ 

+ANN^NZ^,N,N 
+2N^Z^^,NN,N + 2NN^NZ,^,,N^ 



-2N^Z.,NZ.,N^ 



'—R^N- —nl„N 
72 180 



^^RUpN)^ . (5.4) 

As noted, the existence of a Chan's form for the coeffi- 
cients was not completely obvious a priori in the curved 
case. The fact that this Chan's form exists suggests that 
perhaps Chan's method could find a suitable extension 
in the case of curved space. 



VI. THE EFFECTIVE ACTION 

After functional integration, the effective action of a 
complex bosonic field is given by — Tr(logi4r). This can 
be related to the heat kernel by 



Trlogi4:= / — Tr(e^'^). 
Jo T 



(6.1) 



Upon restoring r in the expressions, a contribution 
lri,...,r„0 in Bn picks up a factor T-7+i-p-<i/2^ where 2j 
is the mass dimension of the operator O and p — X^ILi 
(e.g., for l2^2RX^X^ in B2). After carrying 

out the integral over t, the integral over z in Iri,...,r„ 
can be traded by a momentum integral. This gives the 
replacement rule for going from the heat kernel to the 
effective action 

^ ;Ir,,...,r^O ^ IP-Z,r^O (6.2) 



(47r)''/2 



Using p.22|) . An in A^-form is easily obtained from i?„ 



Note that the sign of the fourth term of B2 is incorrect in [Ho 



* Let us note that the minimal parts of A2 in X-form and in N- 
form are different: 
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where we have defined 

. r(d/2) 



new tetrad fields indices. This is achieved through a con- 



d''q 

r(fc + d/2) J (2nf 



{q^YNl^-.-N:,-, (6.3) 



with N = {q^ — X)~^. (Note that /^^^...^^^ depends also 
on d.) The contributions to the effective action may be 
ultraviolet and infrared divergent. Dimensional regular- 
ization applies here. These integrals have been computed 
in (5^ using minimal subtraction. 

The replacement rule gives for the effective action ex- 
panded in derivatives 



with 
Wo 
Wi 
W2 



tr log K = 



n=0 



(6.4) 



(^^2,2,2,2 + '^h,i,3,i)-^t^-^i^^'^^'^ 

~^-^3.3X fj,fj,Xiyi, + 2l2_2^2X fiXuZp^u + 2 "^2,2'^M'"^A"^ 



1 2 2 

In oRX I, X I, — — ^ 

12 M M 3 



1 

30' 



180 



rOT^2 
^2 '^/ji/ 



1 

72' 



rO n2 



(6.5) 



VII. SUMMARY 

We have derived, for the first time, expressions valid 
to all orders in X and to four covariant derivatives for 
the diagonal matrix elements and also for the trace of 
the heat kernel operator in a Riemannian curved man- 
ifold. The expressions presented check previously avail- 
able results for the so-call minimal case ,55] and, when 
reexpanded in powers of X, they reproduce the known 
HMDS coefficients to four derivatives. We also extend 
Chan's formula, originally derived for the effective action, 
to include curvature. As in the flat case, the expressions 
in Chan's form are remarkably simple also in the curved 
case. This simplicity suggests a direct calculation of the 
energy-momentum tensor taking a variation of the effec- 
tive action with respect to the metric. Such a calculation 
has not been addressed here. The method of covariant 
symbols allows to consider more general coordinate con- 
nections, including torsion. This would be of interest 
in the derivation of the Lorentz group generators since 
the coordinate connection couples to them. A virtue of 
our formulas is that they are equally simple for scalar 
wavefunctions and for tensor ones, without redefining the 
gauge connection to include the parallel transport of the 



sistent use of the operators 



which are defined so 



that they are multiplicative. Remarkably the formulas 
obtained are independent of the tensor representation of 
the wavefunction even for the coefficients Bn{x). This is 
because the general formulas for the trace cyclic property 
and integration by parts can also be written in a tensor 
representation independent way (see Appendix IB|) . 
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APPENDIX A: COVARIANT SYMBOLS 

Let us give some details on the calculation of An using 
covariant symbols. Fuller details on the use of covariant 
symbols with derivative expansions can be found in (56j 
for flat space-time and [46] for curved space-time. Using 
the expansions (|3.5p in p.4p . reexpanding the result and 
keeping terms with at most two covariant derivatives, 
gives 



{x\e^\x) 



dz 1 



d'^p 



N - NXa d"N + NXa d^NXa d^N 



-0(V^ 



(Al) 



with = (z—pjj^ — X)^ . (Let us warn that we are using 
a purely imaginary p^, to avoid the proliferation of i's in 
the formulas.) The derivatives with respect to are 
then carried out (using df^N = 2p^N'^). After that, the 
shift z ^ z + allows to isolate the p^ dependence in 

integrals of the type / d'^pe^i^p^^ • • -p^^ which are easily 
evaluated. These steps produce 



{x\e^\x) 



1 



dz 



(4^)d/2 2m 



N + NX^^N^ - 2NX^^N^ 



+2NX^NXf,N'' ~ 
-2NX^,N^Xf,N^ 



ANX^.NX^N" 
3 



+o(v^ 



(A2) 
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with N = (z — X) ^. This expression immediately trans- 
lates into those in p.6p and p.7p . 



APPENDIX B: TRACE AND INTEGRATION BY 
PARTS 

In this work (and in particular in this appendix), all 
formulas hold for arbitrary tensor representations of the 
wavefunctions ip{x) on which the Klein-Gordon operator 
acts. The space of world tensors of rank r is spanned by 
^al ' ' ' ^al where e^(x) is a local basis of the space-time 
tangent space. 

The trace tr in (|2.18p refers to gauge indices and to 
world indices, tr — trgaugotrworid • We need to consider 
only multiplicative operators O (that is, 0-0 at x does 
not depend on ij) dX x' ^ x). All operators in An and 
Bn are multiplicative. In this case, the trace on world 
indices in the representation of tensors of rank r will be 

t^^orM{d)^el\---eX{de>:,l---eZ) (Bl) 

where is the dual basis, e°e{^ = 6ab- Note that 
trworid(C) may only be non vanishing when O is itself 
a world scalar, that is, it maps tensors of rank r to ten- 
sors of rank r. As illustration, in the space of tensors of 
rank r, an easy calculation yields 

tr {^^Zl^j = rf'-trgaugc 

-rd^-^n,^R%^p, (B2) 

Ug being the dimension of the gauge representation. 

The use of the trace cyclic property and integration 
by parts requires some care. These properties work as 



usual when the operators involved are (i) multiplicative 
and (ii) they do not act on world indices. (For a scalar 
operator O, this means that the component /i of O ap- 
plied to ea{x) does not depend on e'^{x), for v ^ ^.) 
Instances of such operators are X, X^ and Vlf^i,. In 
the notation of [1^, these are the operators in the class 
C{Y_,Z_). On the other hand, modifications occur in 
the trace cyclic property and integration by parts when 
multiplicative operators acting on world indices are in- 
volved (class C(V)). An instance of this is Z^j/, since 
{Zfj^uea)^ = Rf^uXaCa + ^f,i^e^- (However, [V^,^] and 
[Z^^,v4] G C(V,Z) provided A e C(V,Z).) 

To write down the correct relations, let us define 
Z^f-iin the curvature parts of Z^^...^^, that is, ob- 
tained by dropping in the covariant derivative. In 
particular, 

Z^ll, = Z^jj + . (B3) 

Furthermore, let us introduce the shorthand notation 

(O) J d'^x^triO). (B4) 

The two following useful properties are easily estab- 
lished 

(AZj!,) = 0, AeCiV,Z) 
{AZ^t^u) = {-\Raac.t..A), AeC(V,Z) (B5) 

which hold for arbitrary multiplicative operators A not 
acting on world indices. 

Using these properties, one can prove the following re- 
lations for arbitrary operators A,B,..., in C(V, Z): 



([V^,A]i3) = {-A[V^,B]), (B6) 

([Vq, ^jBZpjyC) — — A{BaZf^iyC + BZafii/C + BZf^iyCa + Raafii/BCa + -RcrcratJ.uBC)) , (B7) 

(ABZ^.C) - {BZ^,CA-BC[Zi;,„A]), (B8) 

(ABZ^^^C) = {BZ^^,CA-BC[Z^^^,A]), (B9) 

{ABZ^^CZ^fiD) = {BZ^^CZ^pDA- BZ^,CD[Z^p,A]~ BCZ^fiD[Z^^,A]+ BCD[Z^p,[Z^^,A]]) , (BIO) 

{ZfiuAZapB) — (^AZajsBZf^i, — Raf)^j,aABZcy — RafSvaABZ^a) . (Bll) 



The first two identities refer to integration by parts while relation (|4.2p is a consequence of (jBSp . 
the other allow to apply the trace cyclic property. The 
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